Shear-thickening and entropy-driven reentrance 
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We discuss a generic mechanism for shear-thickening analogous to entropy-driven phase reen- 
trance. We implement it in the context of non-relaxational mean-field glassy systems: although 
very simple, the microscopic models we study present a dynamical phase diagram with second and 
first order stirring-induced jamming transitions leading to intermittency, metastability and phase 
coexistence as seen in some experiments. The jammed state is fragile with respect to change in 
the stirring direction. Our approach provides a direct derivation of a Mode- Coupling theory of 
shear-thickening. 
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When liquids are subjected to strong stirring, in gen- 
eral their viscosity decreases, a phenomenon known as 
'shear thinning'. The opposite ('shear thickening') be- 
havior when stirring leads to increased jamming is rather 
exceptional and intriguing That shear-thinning 

should be generic is quite easy to understand by con- 
sidering the structure of any system having a large vis- 
cosity and long relaxation time scales. From the point of 
view of the phase-space energy landscape, the long relax- 
ation times are the consequence of directions that are ei- 
ther almost flat, or contain barriers that can only just be 
crossed. When the system is stirred by non-conservative 
forces, displacements are easily induced along these di- 
rections, and this has the effect of speeding the relax- 
ation. From the real space point of view, slow relaxations 
are linked to extended dynamically correlated spatial re- 
gions [2(, and stirring tends to break off these domains, 
thus making the system more fluid, ft is thus no sur- 
prise that just about any system (and any model) will 
naturally exhibit shear-thinning. 

For shear-thickening instead, several explanations have 
been attempted, and it is at present not clear whether a 
universal one will apply for all possible systems. There 
are indications that jamming in particulate suspensions 
is related to increased disorder [3(, and in some cases 
to the formation of clusters of particles in lubrication 
contact 0, ■ To the extent that at large densities the 
strongly jammed state has the appearance of an amor- 
phous, glassy solid, shear-thickening may be thought of 
as a consequence of an underlying glass transition in- 
duced by stirring HJ. For such collective behavior one 
can attempt a theory with less focus on the details but 
founded on notions that are thought to be generic of 
glasses: this has suggested, for example, the casting of 
the problem in a Mode Coupling format [6j. 

In this Letter we attempt a microscopic setting in 
which glassiness and shear-thickening emerge naturally 
and are simultaneously understood. The basic idea is to 
exploit the analogy between entropy-driven transitions 
in which systems freeze upon heating and those in which 
they jam under the action of stirring. To obtain con- 
crete results, we discuss this idea within the 'Random 



First Order' scenario for the glass transition, although 
it is not restricted to it. In its simplest version, the 
Random First Order scenario applies to models of fully 
connected degrees of freedom and a complicated set of 
interactions with or without quenched disorder 0, @ ■ ^ 
allows for a unified description of the (fragile) glass tran- 
sition, both from the dynamical and from phase-space 
landscape points of view. The approach contains as a 
special, high temperature case the Mode Coupling The- 
ory || , while in the low temperature regime it provides a 
theory for aging 8] . As is characteristic of this approach, 
it has the satisfactory feature that many different aspects 
of the collective behavior follow without further assump- 
tions, and the weakness that spatial features are not (for 
the moment) fully incorporated. 

Introducing stirring. Shear-thinning appears nat- 
urally if one considers the action of 'stirring' terms ca- 
pable of generating permanent currents, i.e. forces that 
do not derive from a (global, time-independent) poten- 
tial. A useful, though approximate, way to see the effect 
of random stirring is the following: consider a system 
with coordinates X{ evolving according to some form of 
dynamics (Langevin, Monte Carlo, molecular dynamics) 
in contact with a heat bath at temperature T and un- 
der the action of a potential and of stirring forces ff " 
acting on the z th degree of freedom. Stirring forces are 
by definition nonconservative, suppose (although this is 



inessential) they are linear: ff" 



If we make 



the simplifying though rather crude assumption that the 
J*!- are long range, randomly distributed (so J|J? is asym- 
metric) and uncorrelated, one can easily show [ToL fll| 
that on average /| tir = p;(i) where pi{t) are Gaussian 
noises with correlations (pi(t)pj(t')) = SijC(t,t') where 
C(t,t') = J2k( x k(t) x k{t')) /N is the two-time autocor- 
relation function. The stirring thus provides a random 
noise unmatched by a friction term: this can be seen 
as a coupling to an infinite temperature (self-consistent) 
bath [12j| . Just like in any stirring situation, if the sys- 
tem for some reason does not flow, the noise pi becomes 
time-independent and hence does no work. 

The generic situation when the (e.g. Monte Carlo) dy- 
namics is perturbed by a nonconservative force is that the 
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structural a-relaxation time becomes shorter - a shear- 
thinning effect - and that a two-temperature regime 
emerges 0] even in the supercooled liquid pha se. This 
scenario has been tested in realistic systems and the 
agreement is impressive. In contrast to the case of shear- 
thinning, there has as yet been no way to introduce or 
understand shear-thickening in these terms, and a phe- 
nomenological construction with a Mode-Coupling flavor 
has to be introduced in a somewhat ad hoc manner, with 
no underlying microscopic model 0|. 

A phase reentrance mechanism. Let us review 
briefly a microscopic mechanism [l^ for freezing induced 
by heating (inverse freezing) . Suppose one has an ensem- 
ble of molecules (e.g. polymers) that have a low temper- 
ature ('folded') state in which they are mutually weakly 
interacting, and a higher temperature ('unfolded') state 
which is favored entropically and in which they interact 
strongly with each other. As temperature is increased, 
each polymer unfolds and reaches out to the other poly- 
mers, the resulting entangling thus may lead to a glass 
transition. A further increase of temperature will even- 
tually lead back to a liquid phase. 

In order to obtain a minimal model of a liquid that 
upon heating is driven by entropy into a glass, one can 
consider spins taking values 0, ±1, and a Hamiltonian 
consisting of a term oc sf favouring the 'folded' con- 
figurations Si — 0, and an interaction term J^ij JijSiSj 
that is active when the spins are in the 'unfolded' states 
Si = ±1: 



(1) 



The entropic favouring of the Si = ±1 configurations 
is enhanced by making these states r-fold degenerate. 
Schupper and Shnerb chose the interactions Jij from a 
fully-connected Gaussian distribution, thus obtaining a 
spin- glass-like phase If instead one wishes to model 
a structural (fragile) glass behavior, one may choose in- 
teractions as in either the random orthogonal model ^3 > 
or to consider a p-spin interaction model with spin-1 vari- 
ables like in Ref. [lj|- We have studied the equilibrium 
phase diagram of the former model in detail and found, 
for large enough r, a reentrant behaviour in both the 
dynamic and static glass transition line. The structural 
glass transition in this model can be either thermody- 
namically first or second order (i.e. with or without la- 
tent heat) , depending on the value of D ^| . 

Shear-thickening models. As mentioned above, 
stirring is somewhat analogous (and in the example above 
exactly equivalent) to coupling to a high temperature 
bath. One can thus imagine that in a problem with phase 
reentrance, stirring might induce a transition from the 
liquid to the glassy phase. This is clear in the folded 
polymer problem described above: taking into account 
the known fact that shearing [2(| , or random stirring |2l| 
can make polymers unfold to the interacting state, this 
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FIG. 1: Energy density vs. time after a system (size TV = 128) 
in the liquid phase at T = 0.03 and D — 3.0 was taken to the 
jammed phase by a stirring force with e = 1.6 and S = 0.2. 



may result in an increase in the viscosity of the polymer 
melt. Shear-thickening is in such cases a form of phase 
reentrance [2*2^ . 

What we have discussed so far suggests that one can 
model shear-thickening by considering a nonconservative 
forcing acting on a reentrant model, for example a force 
field /j Stir acting on the z th spin of QJ: 
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ij °3 > 



(2) 



where — —Jjl and = —Kjl are independent 
Gaussian random variables with zero mean and variance 
1/N. The two components of the force field act indepen- 
dently on the 'folded' and the 'unfolded' configurations 
with stirring strengths s and S respectively |23j |. 

Phase diagram - metastability and coexistence. 
Fig. ^ shows the evolution of the energy of a small sys- 
tem after stirring terms @ were applied to the nonin- 
teracting 'liquid' state. The fraction of ±1 spin becomes 
appreciable and in the energy vs. time plots we observe 
intermittent arrest and flow behaviour with jumps be- 
tween long-lived interacting states. If the temperature 
is low and stirring is not too strong, the system quickly 
falls in a state that is for all practical purposes stable. In- 
creasing the stirring strength the trapping times become 
shorter, and for sufficiently high stirring rates the system 
becomes a normal (non-aging) liquid. We shall not dis- 
cuss in detail the shear-thinning (or rejuvenation) aspect, 
as it has been already extensively discussed in the liter- 
ature 0, Q] . Intermittent situations where the syst em 
jams and unjams as in Fig. have been observed |2 1| . 

The analysis of the temporal evolution of the energy 
and the spin-spin correlation function allows to identify 
different dynamical regimes and construct the dynamical 
phase diagram by varying the several parameters char- 
acterizing the system. Fig. [5] shows a section of such 
a dynamical phase diagram in terms of D and stirring 
strength e. The weakly stirred liquid phase I has a low 
density of interactive sites p ~ 0. The 'jammed' phase II 
is characterized by a fraction of spins in the interactive 
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FIG. 2: A section of the dynamical phase diagram for tem- 
perature T — 0.05, S — 0.0 and r = 6. The transition be- 
tween the phases I and II is first-order below the point A and 
second-order above it. 

±1 state and aging (i.e. the progressive trapping in ever 
deeper interacting states), and the 'liquid' III regime for 
high stirring rates is just the result of shear-thinning of 
the jammed state. For D < 2 and low temperature the 
system is glassy in the absence of stirring. On increas- 
ing D and the stirring force there is 'first order' jamming 
transition (below the point A in Fig. |2J with hysteresis 
in e: along this curve the liquid and jammed phases co- 
exist. Above the point A the transition from phase I to 
phase II is a 'second order' jamming transition without 
hysteresis, but with a regime in which the system forms 
under stress an aging glass. For much larger values of D 
there is no jamming for any stirring, but there is contin- 
uous shear-thickening when the T — e trajectory followed 
passes near a transition line in phase space. 

In these mean-field models the dissipated power scales 
as £ 2 /tv,, where r a is the a-relaxation time of the sys- 
tem [l3|. Comparing to a standard shear flow this sug- 
gests that the amplitude of the driving force, e, plays 
the role of a stress, a, while e/r a is analogous to a shear 
rate, 7. We may thus obtain the standard a vs. 7 flow 
curves by increasing the stirring rate at constant D: they 
turn out to be strikingly similar to those of Ref. . In 
the main frame of Fig. [3] we show an example of such 
flow curves corresponding to the 'full jamming' scenario 
of Ref. [2^|. In this case, one observes an interval of 
stress, in Fig. |3| between 1.65 and 2.2, within which the 
flow rate vanishes, even if the system is ergodic at rest. 
The relaxation time, r a , was estimated as the time in- 
tegral of the normalized spin-spin correlation function. 
Examples of correlation curves are shown in the inset of 
Fig. |21 Shear-thickening behaviour (a slower decay of C) 
is observed when a increases from 1.5 to 2.5, while shear- 
thinning (a faster decay of C) appears for higher stress 
(for a increasing from 2.5 to 4 in the inset of Fig.[3J). 



FIG. 3: Flow curve along the constant D = 4 line for T — .05. 
a = e and 7 = a/r a are the analogue of stress and shear rate 
for a mean-field system. Note how the intersection of the line 
D — 4 with the jammed phase in Fig.^is reflected here. The 
inset shows the corresponding spin-spin correlation function, 
C, vs time t, for increasing stress (system size N — 500). 

In Ref. |26| , an experiment is described in which a con- 
centrated suspension of non-Brownian particles is driven 
by stirring from the liquid to a metastable jammed phase. 
It would be interesting to see experimentally whether 
the opposite situation, when the liquid is the metastable 
phase, may occur. 

Chain fragility and aging of the jammed phase. 
One of the properties of materials that are jammed by 
stirring that we may wish to test in this model is the 
fragility with respect to incremental stresses in a different 
direction [27| • If in a system in phase II we change the re- 
alization of stirring forces from J? a to J? B cos 9 + J-f sin 9, 
and similarly for K^, even for small 9 we find that the 
system responds by rearranging its configuration, the 
faster the larger the value of 9 (9 € [0,7r/2]), see Fig. 

Conclusions. In this Letter we discussed a connec- 
tion between the mechanisms of entropy-driven phase 
reentrance and shear-thickening. This relation may ex- 
ist in some cases just in principle, as the temperatures 
or chemical potentials needed to actually affect substan- 
tially the particles may be in practical situations ex- 
tremely high. The present models are clearly schematic, 
but not much has been put into them and yet we see the 
elementary constituents self-organize to produce stirring- 
induced jamming with aging and intermittency, non- 
Newtonian rheological behaviour like shear thinning and 
thickening, metastability and chain fragility. 

Let us finally mention that one can also construct a 
reentrant continuous model with p-spin interactions H p 
following the same idea outlined above [2^]. From this 
model one can immediately obtain Mode-Coupling with 
reentrance — and also shear thickening by adding stir- 
ring forces of the form J2J to the Langevin dynamics 
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FIG. 4: Chain fragility: mutual correlation, Q(t), of two iden- 
tical jammed systems, in one of which the stirring direction 
has been changed at time t = 0. 9 — tv/2 corresponds to 
the case in which the two systems have independent random 
stirring forces, while in the case 6 = 0.2 the random stirring 
forces are correlated. 

ii = —d Xi H p + rji. Another interesting approach is to 
introduce stirring terms in the Mode-Coupling models of 
colloids with short-range attractive potentials 29] . This 
would allow to investigate the intriguing perspective that 
stirring can drive the liquid-liquid, the liquid-glass and 
also the glass-glass transition in such systems. Work 
along these lines is in progress [28| . 
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